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Constraining white-dwarf kicks in globular clusters 
II. Observational Significance 
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ABSTRACT 

If the winds of an asymptotic-giant-branch stars are sufficiently strong are slightly 
asymmetric, they can alter the star's trajectory through a globular cluster; therefore, 
if these winds are asymmetric, one would expect young white dwarfs to be less radially 
concentrated than either their progenitors or older white dwarfs in globular clusters. 
This latter effect has recently been observed. Additionally the young white dwarfs 
should have larger typical velocities than their progenitors. After phase mixing this 
latter effect is vastly diminished relative to the changes in the spatial distribution 
of young white dwarfs with kicks, so it is more difficult to detect than the change 
in the spatial distribution. The most powerful kinematic signature is the change in 
the eccentricity of the orbits that is revealed through the distribution of the position 
angles of proper motion. 
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1 INTRODUCTION 



ISpruitJ l|f 9981 ) argued that their observed rotation rates of 
white dwarfs may be evidence that they receive mild kicks 
from by asymmetric winds toward th e end of their time 
on th e asymptotic gian t bran ch (AGB) (Vassilia dis fc Wood! 
ll993D . lFellhauer et al.l (|2003h argued that these mild kicks 
could accou nt for a putative dearth of white dwar fs in open 
clusters (e.g. lWeidemarm1ll977l : lKalirai et al.ll200ll ). The dy- 
namics of open clusters provide limited evidence because one 
expects the white dwarfs simply to leave the cluster. One 
can probe the dynamics of AGB winds directly through ob- 
servations of masers, but such attempts are dogged by the 
variability of the star itself which makes it difficult to con- 
strain any relative motion between the centre of mass of the 
wind and that of the star. The globular clusters provide a 
unique probe of white dwarf kicks because both the mag- 
nitude of the expected kick and the velocity dispersion of 
the giants within the cluster are on the order of kilometers 
per second — one would expect a significant signal. On the 
other hand the escape velocity of the cluster may be sev- 
eral times larger so most of the white dwarfs remain in the 
cluster to measure after the kick. The expected signature of 
white dwarf kicks has been observed in M4 and NGC 6397 
|Davis et al.ll2006l ). 

This letter will examine how size of the kick, the initial 
distribution of stars and the size of the observational sam- 
ple affects the expected statistical significance that white 



dwarf kicks have on observations of globular clusters through 
Monte Carlo calculations of the phase-space distribution 
function of stars in the cluster. 



2 CALCULATIONS 



iHevll (|2007h gives the details of the numerical calculations, 
but a quick summary is useful. Clusters of stars can typi- 
cally be modelled with a lowered isothermal profile (or King 
model) (|Michielll 9631 : iKindlf 9661 ) 
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, \P is the gravitational potential, a is the 



characteristic velocity dispersion of the cluster as a whole 
and pi is a characteristic density. 

With time the kinetic energy within the cluster ap- 
proaches equipar tition betwee n the various stars such that 
rrnoi = mjaj (jSpitzerl If 9871 ). The progenitors of young 
white dwarfs will be the most massive main-sequence stars 
in a cluster at the time, so they will typically have oto < 
where a is the mean velocity dispersion of the cluster; fur- 
thermore, because these progenitors will only have a small 
fraction of the mass of the cluster, they can be considered 
as massless tracers; their phase space density will be given 
by Eq. {T} with a — » <jto- 
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During their time on the AGB, the stars may lose 
a lar ge fraction of their mass suddenly and asymmetri- 
cally (fvassiliadis fc Woodlll993l : ISpruitll998l : iFellhauer et~aH 
2003) and receive an impulsive kick. Because the distribu- 
tion no longer depends on constants of the motion alone, the 
distribution function itself will depend on time, so a Monte 
Carlo realization of the kicked distribution function is help- 
ful to make further progress. Specifically, for a variety of 
values of *(0), ctto and cr^, one hundred thousand stars are 
drawn from Eq. and given a three-dimensional velocity 
kick drawn from a Gaussian of width of at- Finally, each 
star is evolved forward along its orbit to a random phase. 

The result of the Monte Carlo realization is a list of 
positions and velocities just before the kick, just after the 
kick, and at a random time later. For each star these posi- 
tions and velocities are projected along a random direction 
and into the plane perpendicular to that direction, yielding 
a distribution of velocities and positions parallel and per- 
pendicular to the line of sight to the globular cluster. In 
principle all three components of the velocity and position 
are measurable. A useful measurement of the distance of 
a faint white dwarf along the line of sight requires micro- 
arcsecond astrometry; measurements of the velocity along 
the line of sight require multiplexing spectroscopy again to 
faint magnitudes. Although these may be feasible in the fu- 
ture, the focus will be on the quantities perpendicular to the 
line of sight. Finally, because the globular cluster is spheri- 
cally symmetric the position angle of a star relative to the 
centre of the cluster is uniformly distributed. However, the 
position angle of the proper motion relative to the projected 
direction toward the centre of the cluster does carry useful 
information about the properties of the stellar orbits. 



3 RESULTS 

This leaves five quantities of interest: the distance and speed 
of a star relative the centre of the cluster in the plane of the 
sky, the angle between the proper motion of the star and the 
direction toward the centre of the cluster and of secondary 
interest, the distance and speed of a star relative the centre 
of the cluster along the line of sight. Furthermore, because 
the line of sight to the observer is arbitrary the distribu- 
tions of these latter two quantities are simply related to the 
distributions of the astrometric quantities so the results will 
focus on the three astrometric measurements. 



3.1 Distributions 

The kicks that the white dwarfs receive as they form can 
dramatically affect their phase-space distribution. The up- 
per panel of Fig. [1] shows that a kick of the same mag- 
nitude as the velocity dispersion of the population signifi- 
cantly puffs up the projected radial distribution of the young 
white dwarfs. Additionally, the effect is similar in magni- 
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for the three- 



This agrees with the results of [ 
dimensional radial distribution. 

The more subtle effect is the change in the distribution 
of proper motions, both in magnitude and direction. The 
lower panel of Fig. [T] shows that the change in the distribu- 
tions of the magnitude proper motions is smaller especially 
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Figure 1. The cumulative position and velocity distributions in 
the plane of the sky for uxo = = 0.5<r (left curves) and 
trxo = u k = Q.8<7 (right curves) with VE'(O) = 6a 2 . The solid 
curves follow the initial distribution, and the short-dashed curves 
follow the distribution immediately after the kick (these two co- 
incide for the radial distribution). The long-dashed curves follow 
the distribution after phase mixing. 



for the loosely bound progenitors (ctto = 0.8ct). Immedi- 
ately after the kick the velocity distributions of the young 
white dwarfs spread out dramatically. However, after averag- 
ing over the phase of the new orbits, the effect is diminished. 
Although this may initially be surprising, it results from the 
negative specific heat of gravitationally bound systems. For 
example, if the orbit were initially circular and Keplerian, 
the kick would serve to lift the apocentre of the orbit; there- 
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Figure 2. The cumulative distribution over uniform of the angle 
between the proper motion and the projected direction toward 
the centre of the cluster for (Tto = 0.5<r (solid curves) and 0.8ct 
(dashed curves) with *f(0) = 6<r 2 after phase mixing. From bot- 
tom to top the curves trace a k = 0.5, 1 and 1.5 oto- The distribu- 
tions before the kick and before phase mixing are consistent with 
a uniform distribution. All the distributions for cr k = are also 
uniform. In all cases a uniform distribution has been subtracted 
from the cumulative distribution. 



fore, on average the velocity along the orbit would decrease. 
Of course, the orbits are not initially circular and thew po- 
tential is not Keplerian so the net result is to diminish the 
effect (and not reverse it). The reduction is more dramatic 
for the loosely bound stars because a modest increase in the 
velocity can greatly increase the size of the orbit and possi- 
bly unbind the star from the cluster; therefore, the stars that 
remain bound to the cluster will necessarily have lower ve- 
locities — the loosely bound stars are more strongly affected 
by evaporative cooling. 

The orbits of the stars change in a second way. The or- 
bits after the kick and orbit averaging are typically more ec- 
centric that before. The distribution of velocities according 
to Eq. fT} is spherically symmetric. Furthermore, the kicks 
are distributed symmetrically; therefore, the distribution of 
the angle between the proper motions and the projection 
direction to the centre of the cluster is uniform. However, 
after the kicked stars have mixed in phase, it is more likely 
to find stars moving within forty-five degrees of the radial 
direction than within forty-five degrees of the tangential di- 
rection - the orbits are typically more eccentric than before. 
Fig. [2] shows the cumulative distribution of this angle (0) 
for various sizes of kick with a uniform distribution sub- 
tracted away. First, the effect weakens with the size of the 
kick and vanishes for no kick. Second, for large kicks the 
effect is weaker for the more loosely bound stars. Again this 
is because the loosely bound stars are more likely to get 
eccentricities greater than unity and leave the system. 



3.2 Detectability 

Although these differences in the distributions are apparent 
when many thousands of stars are use as probes, in real- 
ity one typically has only several tens of stars to charac- 
terise the distribution of young white dwarfs (the kicked 
stars) and possibly a factor of ten more stars to estimate 
the distribution of the progenitors. To estimate the typi- 
cal sensitivity of measurements to the effects of white dwarf 
kicks, smaller samples were generated with various numbers 
of kicked stars and a factor of ten more progenitors for com- 
parison. A total of one thousand samples were generated in 
each group with <tto = 0.5a and 0.8a and au = ctto- For 
each pair of samples the likelihood that they were drawn 
from the same underlying distribution was calculated using 
the Kolmogorov-Smirnov test and the Wilcoxon rank-sum 
test, the final result is the median likelihood over the one 
thousand realisations. For the sake of clarity and definite- 
ness observational errors have not been included in the sam- 
ples generated - such errors would reduce the significance of 
the detections. 

The statistical likelihoods depicted in Fig. [3] of detect- 
ing the change in the distributions support what is apparent 
from the distributions themselves. Finding a change in the 
radial distribution is much more straightforward than find- 
ing changes in the velocity distribution (even with perfect 
data). Furthermore, for the loosely bound progenitors (lower 
panel of Fig. [3j , if it is impossible to distinguish the proper 
motion distribution of the young white dwarfs with kicks 
from their progenitors with fewer than one white dwarfs in 
the sample (and 1,000 turn-off stars). The key to finding 
these changes is building up a sample of many stars. For a 
kick comparable to the initial velocity dispersion, at least 
one hundred stars are required even in the best case for 
tightly bound progenitor stars (upper panel of Fig. [3]). 

3.3 Predictions 

iDavis et al.l (|2006l ) reported a significant difference between 
the distribution of young white dwarfs and turn-off stars 
and between young white dwarfs and older white dwarfs. 
Furthermore, they estimated that the distribution of young 
white dwarfs was most similar to that of main-sequence stars 
whose mass was half that of stars near the turn-off. iHevll 
(2007) found that such a large change in the stellar distri- 
bution would require a kick with at > uto- If one starts 
with a to = 0.5a, a kick of 1.82<tto is required to make the 
final distribution of projected radii look like a population 
with dispersion of 0.7a; that is one that corresponds to one- 
half the mass of the progenitors. The difference between the 
initial and final distribution of radii is detectable with a me- 
dian likelihood of the null hypothesis of less than one part 
in one thousand with only twenty young white dwarfs in the 
sampl e, comparable to the probability quoted bv lDavis et al.l 
(|2006h . 

With about forty young white dwarfs, the difference 
between the initial and final proper motion distributions in 
Fig. [4] is detectable with the likelihood of the null hypoth- 
esis of less than two percent. However, the systematics of 
detecting an increase in the typical proper motion of the 
faint young white dwarfs relative to the much brighter main 
sequence stars is problematic because the expected errors in 
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Figure 3. The median significance of a detection of a differ- 
ence between the distribution of young white dwarfs and main- 
sequence stars near the turn-off. For the upper panel the initial 
distribution had cr-po = 0.5<x and a oj. = <t t q. For the lower 
panel the initial distribution had ctto = 0-8cr and a <r fe = ctxo- 
The cluster as a whole has ^(O) = 6ct 2 . The dashed curves follow 
the median likelihood that the two realisations differ by chance 
(the probability of the null hypothesis) over 1,000 trials using the 
Kolmogorov-Smirnov test and the solid curves give the results us- 
ing the Wilcoxon test. The number of main-sequence stars in the 
sample is assumed to be ten times the number of white dwarfs. 
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Figure 4. The proper motion distr ibution for a whi te-dwarf 
kick large enough to account for the Da vis et al. (2006) results 
(cr fe = 1.84<tto> °TO = 0.5cr and *&(0) = 6<r 2 ). In the main figure 
(similar to the bottom panel of Fig. [TJ, the solid curves follow the 
initial distribution, the short-dashed curves follow the distribu- 
tion immediately after the kick, and the long-dashed curves fol- 
low the distribution after phase mixing. The dot-dashed curves 
follow the initial proper-motion distribution for a dispersion of 
0.7a. The inset shows the distribution of position angles of the 
proper motion (similar to Fig- [2J. The solid curve gives the dis- 
tribution if one neglects measurement errors and the dashed lines 
assumes a measurement error in the proper motion of a in each 
direction. From Tab.[T]this translates to an accuracy of about one 
milliarcsecond per year for many nearby clusters. 



the astrometry increase as the stars become fainter which is 
just the effect predicted for white-dwarf kicks. 

On the other hand, an exploration of the direction of 
the proper motion should not be affected by the system- 
atics of errors in the astrometry because the errors in the 
proper motion should be symmetric and even if they are not 
symmetric it is unlikely that they will correlate with the di- 
rection toward the centre of the cluster. For the parameters 
used here, about sixty percent of the stars have proper mo- 
tions within forty-five degrees of the radial direction, and 
about one hundred young white dwarfs would be required 
to detect the asymmetry with the median likelihood of the 
null hypothesis (Wilcoxon rank-sum text) of less than two 
percent. Although more white dwarfs are required for this 
test, it is less plagued by systematics, so in practice it should 
be more powerful. 

In one includes a measurement error in the proper mo- 
tion of a about and assumes that the error is isotropic, 
this necessarily weakens the asymmetry. On the other hand, 
these errors even if poorly estimated cannot create a sig- 
nal where there is none, in contrast with measurements 
of the magnitude of the proper motion. From Tab. \T\ this 
error translates to about one milliarcsecond per year for 
nearby globular clusters, similar to the precision quoted 
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Table 1. Kinetic parameters and distances for a few nearby globular clusters whose radial stellar distribution could be probed with JWST 
(Riche r priv. comm.). The concentration or the ratio of the core radius to the tidal radius (c = log 1n rt/r c ) and the dista nce are from lHarrij 
119961 1, a dash denotes that the core has collapsed. The values of ^(O)/^ are fromlMcLaughli n &c van der Marell ||2005|1 with the exception 
of NGC 6837 where it was inferre d fro m the concentration . Th e value of a c was obtained by as suming that M/Mxo = 0.3/0.8 = 0.375. 
Refs: fal lMcLaughlin et alj j2006t) ; fbl lReiins et al.l d2006t) ; (c) iPeterson fc Lathaml dl986T) ; (dl lPrvor fe Mevlanl dl993r i. 



Cluster 


<r los [km/s] 


a [km/s] 


c 


*(0)/o2 


d [kpc] 


cr/d [milliarcseconds/yr] 


NGC 104 (47 Tuc) 


11.6±0.8( a ) 


19 


2.03 


8.60±0.10 


4.5 


0.89 


NGC 5139 (Omega Cen) 


15(b) 


25 


1.61 


6.20±0.20 


5.3 


1.00 


NGC 6121 (M4) 


3.9±0.7< c ) 


6.4 


1.59 


7.40±0.10 


2.2 


0.61 


NGC 6397 


3.5±0.2< d ) 


5.7 






2.3 


0.52 


NGC 6656 (M22) 


8.5±1.9< c ) 


14 


1.31 


6.50±0.20 


3.2 


0.92 


NGC 6752 


4.5±0.5< d > 


7.3 






1.0 


0.38 


NGC 6809 (M55) 


4.2±0.5< d) 


6.8 


0.76 


4.50±0.10 


5.3 


0.27 


NGC 6838 (M71) 


2.8±0.6( c ) 


4.6 


1.15 


5.4 


4.0 


0.18 



bv lRicher et a"D l|2007h . Including these realistic errors one 
would expect about 52.5 percent of the stars to have proper 
motions within forty-five degrees of the radial direction, so 
about 1,600 stars would be required to detect this asymme- 
try at the two-sigma level if the astrometric errors in current 
instruments are included. 



4 DISCUSSION 

If white dwarfs do indeed receive kicks as they form, both 
their radial and velocity distributions will appear different 
from their progenitors and the older white dwarfs whose 
distribution has relaxed. The easier of these effects to mea- 
sure is the chang e in the radial distribution of the stars 
l|Davis et al.l 120(361 ) : however, there are two other lines of 
evidence that might be present in astrometric measure- 
ments of globular clusters. First, white dwarfs typically have 
larger proper motions than their progenitors and older white 
dwarfs. Unfortunately, the errors in the astrometry of young 
white dwarfs are much larger than their much brighter pro- 
genitors on the main sequence, so unless this systematic is 
handled carefully, this signal is difficult to unravel. Second, 
the white dwarfs after the kick are more likely to lie on ra- 
dial orbits than their progenitors. Although measurement 
errors will weaken observations of this effect, they cannot 
mimic an asymmetry in the orbit parameters, so this sig- 
nal although weaker is less prone to systematics and may 
eventually provide an independent line of evidence for white 
dwarf kicks. 

Furthermore, these calculations indicate that we are on 
the verge of being able to probe the dynamics of globular 
clusters in an entirely new way through the direct measure- 
ment of the proper motions of the stars. Such measurements 
will provide unique information on the structure of clusters 
as well as stellar evolution. 
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